Reconstructing the Cosmic Equation of State from Supernova distances 
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Observations of high-redshift supernovae indicate that the universe is accelerating. Here we present 
a model-independent method for estimating the form of the potential V((j)) of the scalar field driving 
this acceleration, and the associated equation of state w^. Our method is based on a versatile ana- 
lytical form for the luminosity distance Dl, optimized to fit observed distances to distant supernovae 
and differentiated to yield V{cj)) and w^. Our results favor w^ ~ —1 at the present epoch, steadily 
increasing with redshift. A cosmological constant is consistent with our results. 
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The observed relation between luminosity distance and 
redshift for extragalactic Type la Supernovae (SNe) ap- 
pears to favor an accelerating Universe, where almost 
two-thirds of the critical energy density may be in the 
form of a component with negative pressure ||l|-^. Al- 
though this is consistent with J^m < 1 and a cosmolog- 
ical constant A > {e.g. ||]), at the theoretical level a 
constant A runs into serious difhculties, since the present 
value of A is '--^10^^'^ times smaller than predicted by most 
particle physics models 0. 

However, neither the present data nor the theoreti- 
cal models require A to be exactly constant. To explore 
the possibility that the A-like term {e.g. quintessence) is 
time-dependent, we use a model for it that mimics the 
simplest variant of the inflationary scenario of the early 
Universe. A variable A-terni is described in terms of an 
effective scalar field (referred to here as the A-field) with 
some self- interaction V{(l)), which is minimally coupled 
to the gravitational field and has little or no coupling to 
other known physical fields. In analogy to the inflation- 
ary scenario, more fundamental theories like supergravity 
or the M-theory can provide a number of possible can- 
didates for the A-field but do not uniquely predict its 
potential V{4)). On the other hand, it is remarkable that 
V{4>) may be directly reconstructed from present-day cos- 
mological observations. 

The aim of the present letter is to go from observations 
to theory, i.e. from Dl{z) to V{4)\ following the pre- 
scription outlined by Starobinsky Q (see also p). This 
is the first attempt at reconstructing V{<j)) from real ob- 
servational data without resorting to specific models {e.g. 
cosmological constant, quintessence etc.). 

Since the spatially flat Universe {Q,^ + Q,m = 1) is 
both predicted by the simplest inflationary models and 
agrees well with observational evidence, we will not 
consider spatially curved Friedmann-Robertson- Walker 
(FRW) cosmological models. In a flat FRW cosmology, 
the luminosity distance D^ and the coordinate distance 
r to an object at redshift z are simply related as (c = 1 



here and elsewhere) 



a^r = ao 



dt' 



Dl{z) 



n a{t') l + z' 
This uniquely defines the Hubble parameter 
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Note that this relation is purely kinematic and depends 
neither upon a microscopic model of matter, including a 
A-term, nor on a dynamical theory of gravity. 

For a sample of objects (in this case, extragalactic 
SNe la) for which luminosity distances Dl are mea- 
sured, one can fit an analytical form to D^ as a func- 
tion of z, and then estimate H{z) from (0). If p„j — 
{3Hq /8TrG)i^M{a/ao)^^ is the density of dust-like cold 
dark matter and the usual baryonic matter, then 



H' = gTTG ( p„ 



1 



+ 7;r + V{(b) 



from where it follows that 



H^^47rG{pra+4''')- 



(3) 



(4) 



Eqs. (Q) & (|j) can be rephrased in the following form 
convenient for our current reconstruction exercise, 
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dx 



ZH^x 



dx 
dlnH flux 

1p~ 



^^mx^, 



dx 



(5) 
(6) 



where x = l + z. Thus from the luminosity distance D^, 
both H{z) and dH{z)/dz can be unambiguously calcu- 
lated. This allows us to reconstruct the potential V{z) 
and d(f>/dz if the value of flM is additionally given. In- 
tegrating the latter equation, we can determine 4>{z) (to 
within an additive constant) and, therefore, reconstruct 



the form of V{(j)). Note also that the present Hubble 
constant Hq = H{z = 0) enters in a multiplicative way 
in all expressions. Thus, neither the potential V{(J))/Hq 
nor the cosmic equation of state w^{z) depends upon the 
actual value of _ffo- 
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FIG. 1. The maximum deviation /^Dl/Dl between the 
actual value and that calculated from the ansatz (M) in the 
redshift range z =0-10, as a function of r2M = 1 — Q.^. The 
three curves plotted are for constant values of the equation of 
state parameter (as dehned in Eq. [L3) w^ = — 1 (solid line), 
—2/3 (dotted line) and —1/3 (dashed line). 



A fitting function for D^ '■ We use a rational (in terms of 
y/x) ansatz for the luminosity distance D^, 



X 



2 



X — a^Jx — 1 + a 



iix + ^\fx + 2 — a — /3 — 7 



(7) 



where a, /? and 7 are fitting parameters. This function 
has the following important features: it is valid for a wide 
range of models, and it is exactly equal to the analytical 
form given by (1^) for the two extreme cases: Q^ — 0, 1. 
At these two limits, as 57m ^ 1, a + 7 ^ 1 and (3^1; 
and as 51m "*■ Oj ct, /?, 7 -^ 0. The accuracy of our ansatz 
is illustrated in Fig. ^ 

We choose this form since the value of H{z) obtained 
by differentiating Dl/x, according to (g), has the cor- 
rect asymptotic behavior: H{z)/Ho — > 1 as z ^ 0, and 
H{z)/Ho = nll^(l + z)3/2 for z > 1, where 



flu — 



(3' 



a/3 + 7 



(8) 



This ensures that at high-z, the Universe has gone 
through a matter dominated phase, ft should be noted 
that f2M can be slightly larger than the CDM compo- 
nent Qm since the A-field (or quintessence) can have an 
equation of state mimicking cold matter (dust) at high 
redshifts. For instance, Qm — 1-1 ^m in the quintessence 
model considered by Sahni & Wang [nO[. On the other 
hand, JIm ^ 1-15 ilM, to ensure that there is sufficient 
growth of perturbations during the matter-dominated 
epoch (see, e.g., the relevant discussion in [pj). 



Note that the right hand side of (^) should be non- 
negative for the minimally coupled scalar field model. 
At z = 0, this condition gives 



4/3 -f 27 - a 
2-a 



> 3n 



M, 



(9) 



where the equality sign occurs when the A-term is con- 
stant. The fact that D^ is smaller in a universe with 
time-dependent A-term than it is in a constant-A uni- 
verse leads to a lower limit for the parameter /3. When 
taken together with the fact that /? ^ 1 as Hm -^ 1 
(O^ -^ 0) this leads to the following set of constraints 
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FIG. 2. The distance modulus (m - M) of the SNe la 
relative to an ^m — > Milne Universe (dashed line) , together 
with the best-fit model of our ansatz (H), plotted as the solid 
line. The extreme cases of the {Qm, ^4,)— (0, 1) and (1, 
0) universes are plotted as dotted lines. Also plotted as the 
dot-dashed line is the best fit Perlmutter et al. M model (f^M, 
0^)= (0.28, 0.721 The filled circles are the 54 SNe of the 
"primary fit" of |l[. The high-z SNe of || (not used in this 
analysis) are plotted as open circles. 



The observational data: Till date, about 100 SNe la in 
the redshift range z= 0.1—1 have been discovered, a large 
fraction of which have reliable published data from which 
luminosity distances can be calculated. We use the 54 
SNe la from the preferred "primary fit" ('C in their Ta- 
ble 1) of the Supernova Cosmology Project m, including 
the low- 2: Calan Tololo sample Q as used therein. We 
adopt the quoted redshifts, reducing them to the cosmic 
microwave background frame. 

Maximum likelihood fits: The luminosity distance Dl 
(Mpc) is related to the measured quantity, the corrected 
apparent peak B magnitude ms as ms = Mq + 25 + 
5 logj^Q Dl , where Mq is the absolute peak luminosity of 
the SN. The function to be minimized is 






[yizi) - y{mBi)Y 



y(z) EE 10^»/5Di(z). (11) 



TABLE I. Best-fit parameters 
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0.25 



0.3 
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1 702+"-l'^l 
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-0 2S0+°''*^ 

1.226 ±0.018 

1.03 

0.2 ±0.11 



1 630+"-^"'^ 

604+°°^^ 
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U.ZOO^o Q52 

1.229 ±0.018 
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1 cqq+O.iiUa 

±.ooo_o 575 

U.OiO_o 027 
n oro+O.OSe 
U.ZOZ_o 061 

1.232 ±0.018 

1.03 

0.3 ±0.23 



''(X^> = X?nin/(A^-6), where iV = 54. The fit of four variables 
is subject to the constraint CIm ~ ^m, and the inequality (bl), 
which significantly restricts the permitted region. 

A fourth fitting parameter, k = 2x 10*^"/^(c/i/o), which 
is required in addition to a, /3, 7 in the above mini- 
mization process, includes both Mq and Hq, which can- 
not be measured independent of each other. For in- 
stance, if Mq = —19.5 ± 0.1 and flu = 0.3, the value 
of iJo = 61.3±2.9 km s^^ Mpc^^. Note that k only fea- 
tures in the fit of (R) to the data, and does not play a 
role in the reconstruction of V{(j)). 

To obtain the best fit model, we perform an orthogonal 
chi-square fit, using errors on both the magnitude and 
redshift axes in Ui, subject to the constraints (y), (|lO|) and 
the condition J7m — ^m- The latter condition is used for 
simplicity - our results remain essentially the same even 
if we use the entire permitted range VIm ^ l.lSilM- 

The results shown in Table | and in Figure g are for 
f^M = 0.3. In arriving at the best fit, the two constraints 
in ( [L0[ ) are found to be redundant, which means that only 
two constraints, <^) and f^M = ^m, are actually used. 
Reconstructing the scalar field potential: We show the 
form of the effective potential V{z) reconstructed using 
(||) in Fig. H along with the corresponding plot for V^(0), 
where is calculated by integrating (H). The field (f> is 
determined up to an additive constant (J)q, so we take (f> 
to be zero at the present epoch (2 = 0). 

Our experiments with several realizations of synthetic 
data show that this method works best if we fix the 
value of r^M- Henceforth, all reconstructed quantities are 
shown for Qm = 0.3. 

For a scalar field, the pressure p = —T^ — i(/)^ — V 
and the energy density e = Tq = ^(p^ + V are related by 
the equation of state. 



w^{x) = - 



p {2x/3)dlnH/dx-l 



1 - (Hi/m) nux-' 



(12) 



For the Cosmological constant, w = —1, while 
quintessence models [|l2| generally require — 1 < w < 
for z<2. 

Our reconstruction for w^{z) according to dl2) is plot- 
ted in Fig. 1^. There is some evidence of possible evolu- 
tion in W0 with — l<u'0<— 0.86 preferred at the present 
epoch, and — 1 < W0 ^ —0.66 at z = 0.83, the farthest 
SN in the sample (both at 68% confidence, upper hmits 
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FIG. 3. The eflective potential V{z), and the kinetic en- 
ergy term <j)'^, are shown in units of pa = SHq/SttG. Also 
plotted are the two forms of V{(j>) for this V{z), where the 
errors do not reflect errors in the z-ff) relation. The value of 
(known up to an additive constant) is plotted in units of the 
Planck mass mp. The solid line corresponds to the best-fit 
values of the parameters. In each case, the shaded area cov- 
ers the range of 68% errors, and the dotted lines the range 
of 90% errors. The hatched area represents the unphysical 
region t/)^ <0. 



correspond to —0.80 and —0.46 at 90% confidence respec- 
tively). However, a cosmological constant with w = —1 
is consistent with the data. 




Redshift z 

FIG. 4. The equation of state parameter W4,{z) = P/ p as a 
function of redshift. The solid line corresponds to the best-fit 
values of the parameters. The shaded area covers the range 
of 68% errors, and the dotted lines the range of 90% errors. 
The hatched area represents the region w^ < — 1, which is 
disallowed for a minimally coupled scalar field. 

The errors quoted in this paper are calculated using a 
Monte-Carlo method, where, in a region around the best- 
fit values of the parameters shown in Table 1, random 



points are chosen in parameter space from the probabil- 
ity distribution function given by the x^-function that is 
minimized to yield the best fit. At each value of z in the 
given range, the function in question is evaluated at over 
10^ such points, and the errors enclosing 68% and 90% 
of all the values centered on the median are shown in the 
figures. 

The ages of objects: Our ansatz (0) also provides us with 
a model-independent means of finding the age of the uni- 
verse at a redshift z, 



i(z) = H, 



dz' 



il + z')h{z')'' 



(13) 



where the value of h{z) = H{z)/Hq is determined from 
(0) . Figure S shows the age of the Universe at a given z 
and compares it with the ages of two high redshift galax- 
ies and the quasar B1422-I-231 |jl||. We find that the 
requirement that the Universe be older than any of its 
constituents at a given redshift is consistent with our 
best-fit model, which is a positive feature since a flat 
matter-dominated Universe must have an uncomfortably 
small value of Hq to achieve this. 



the potential of a minimally coupled scalar A-field (or 
quintessence) and its equation of state w^{z). It should 
be noted that the basic equations of this ansatz: (g), (|j), 
( p^ ) & (|3|) are flexible and can be applied to models 
other than those considered in the present paper. For in- 
stance one can venture beyond minimally coupled scalar 
fields by dropping either one or both of the constraints 
(^) & ( p^ ) (this is equivalent to removing the constraint 
PA + PA > on the A-field) . Even with the limited high- 2; 
data currently available, our ansatz gives interesting re- 
sults both for the form of V{(j)) as well as Wcf,{z). As data 
improve, our reconstruction promises to recover 'true' 
model-independent values of T^((/>) and w^{z) with un- 
precedented accuracy, thereby providing us with a deep 
insight into the nature of dark matter driving the accel- 
eration of the universe. 
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FIG. 5. The age of the Universe at a redshift z, given in 
units of Hq^ (left vertical axis) and in Gyr, for the value of 
Ho = 61.3 km s~^ Mpc~^ (right vertical axis). The shaded 
region represents the range of 68% errors, and the dotted 
lines the range of 90% errors. The three high-redshift objects 
for which age-dating has been published llSl are plotted as 
lower limits to the age of the Universe at the corresponding 
redshifts. The dashed curve shows the same relation for an 
(f^M, f2^) = (l,Q) Universe for the same Hq. 

Discussion: In this letter, we have proposed a simple, 
analytical, three parameter ansatz describing the lumi- 
nosity distance as a function of redshift in a flat FRW 
universe. The form of this ansatz is very flexible and 
can be applied to determine Dl either from supernovae 
observations (as we have done) or from other cosmolog- 
ical tests such as lensing, the angular size-redshift rela- 
tion etc. Using the resulting form of Dl we reconstruct 
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